Angular dependence of the melting of the porous vortex matter in Bi 2 Sr 2 CaCu20 
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Vortex matter in Bi2Sr2CaCu2 0s with a low concentration of tilted columnar defects (CDs) 
was studied using magneto-optical measurements and molecular dynamics simulations. It is found 
that while the dynamic properties are significantly affected by tilting the magnetic field away from 
the CDs, the thermodynamic transitions are angle independent. The simulations indicate that 
vortex pancakes remain localized on the CDs even at large tilting angles. This preserves the vortex 
thermodynamics, while vortex pinning is considerably weakened due to kink sliding. 



The dynamic and thermodynamic properties of vortex 
matter in high T c superconductors are strongly affected 
by the introduction of correlated disorder in the form of 
columnar defects (CDs). In the presence of a low con- 
centration of CDs, experimental [l|, [2|, |3[ and theoretical 
results 0, 0, 0, 0, B 0] indicate the presence of two vor- 
tex subsystems: vortices localized on pinning centers and 
interstitial vortices between them. At low temperature, 
a porous vortex solid phase is created. A rigid matrix of 
vortices, strongly pinned on a network of random CDs, is 
formed, while the interstitial vortices create relatively or- 
dered nanocrystals within the pores of the matrix. When 
the temperature is increased, the interstitial vortices melt 
while the rigid matrix remains localized, forming a nano- 
liquid phase in which nanodroplets of liquid are weakly 
confined within the pores of the matrix. When the tem- 
perature is increased further, the matrix delocalizes at 
the derealization transition, and the liquid becomes ho- 
mogeneous. Above this transition the effect of the CDs 
on vortex dynamics becomes indiscernible. 

This picture holds when the angle between the av- 
erage induction B and the CDs is zero. According to 
theoretical models [Io|, EH, for small the localized vor- 
tices on CDs remain completely trapped along the CDs. 
At larger the vortices form a staircase structure with 
segments trapped on different CDs connected by weakly 
pinned kinks. At still larger angles, the vortices follow the 
field direction and are not affected by the correlated na- 
ture of the CDs. In general, one expects the effect of the 
correlated nature of the CDs to decrease as 6 increases, 
and to vanish at = 90°. In more isotropic materials 
such as YBa2Cu307, this picture is supported by var- 
ious experiments, showing an angle dependent behavior 
in both dynamic [13] and thermodynamic [l3[ properties. 
In layered materials such as Bi2Sr2CaCu208 (BSCCO), 
vortices consist of stacks of two-dimensional pancake vor- 
tices (PVs) defined in the superconducting Cu02 layers 
only. Measurements of the equilibrium torque and mag- 
netization in BSCCO Q have shown that the CD occu- 
pation by PVs in the liquid phase is angle independent. 



In contrast, Josephson plasma resonance measurements 
[15[ show that the PV alignment in the liquid is signifi- 
cantly affected by CDs. The corresponding equilibrium 
properties of the solid phase have not been studied. Since 
the melting transition depends on the difference in the 
behavior of the free energies of the solid and liquid, it 
is very interesting to study the angular dependence of 
the melting in the presence of CDs. Moreover, CDs are 
known to shift the first-order melting line upwa rd ll| and 
point disorder was found to shift it downward [16]. One 
might expect, therefore, that tilting B away from the 
CDs would shift the melting line continuously from above 
the pristine melting line to below it, reaching the lowest 
position at = 90°, where CDs act like point defects. 

Here, we study the angular dependence of the dynamic 
and thermodynamic properties of porous vortex matter 
in BSCCO. Our measurements show that the irreversibil- 
ity line shifts to lower T when the field is tilted away from 
the CDs, as expected H, [l7]. However, the thermody- 
namic melting and derealization lines have no angular 
dependence. Our simulations show that even at large an- 
gles the PVs are effectively pinned by the CDs and there- 
fore the thermodynamic features remain unaffected. On 
the other hand, the formation of a kink structure of the 
PV stacks increases their mobility due to kink sliding, 
therefore suppressing the irreversible properties. 

For the measurements, BSCCO crystals (T c « 90 K) 
were irradiated, at GANIL, with 1 GeV Pb ions, through 
stainless steel masks with triangular arrays of 90 jam 
holes. This produces amorphous columnar tracks, with 
density rid equal to the ion dose, only in those parts of 
the crystals situated under the holes. Crystals A and 
B, irradiated at 45° with respect to the c-axis, have 
matching fields = n^o of 20 G and 30 G respec- 
tively. The measurements were performed using differen- 
tial magneto-optical imaging (DMO) [18] with a vector 
magnet that allows application of a DC magnetic field 
in any direction, so that the local magnetic response can 
be imaged over the full range of angles from B || CDs 
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FIG. 1: (color online) (a) The irreversibility line Hf? r (T) of 
crystal A, mapped for five angles between B and the CDs, 
from = to 90°. is the pristine melting line, (b) DMO 
image of part of the sample at T = 79 K and H z = 35 Oe. 
(c) Flux profiles across the irradiated region [dashed line in 
(b)], at 67 K and H z = 22 Oe, for = 0° and 90°. 



((9 = 0) toB±CDs ((9 = 90°). 

Figure 1(a) shows the irreversibility line of the irradi- 
ated regions H^ r (T) on an H z —T plane, where H z is the 
c-axis component of the applied field, for different angles 
of B with respect to the CDs [19]. The curves were ob- 
tained using DMO with field modulation of 1 Oe, upon 
sweeping the temperature. Below H^ r (T), the exter- 
nal field modulation is shielded in the irradiated regions 
due to enhanced pinning [l|. Hence, the dark regions in 
the DMO image (Fig. 1(b)) are the irradiated regions, 
while the bright background is the reversible pristine re- 
gion. Upon increasing T, these dark regions gradually 
fade until they disappear completely at H^ r (T), where 
the screening current drops to an unobservably low level. 

Figure 1(a) shows that upon tilting the field away from 
the CDs (changing 6 from to 90°), the irreversibility line 
shifts to lower temperatures. For example, at H z = 70 
Oe the irreversibility temperature at 90° is 2.5 K lower 
than at = 0. Even at = 90°, H?* r (T) remains slightly 
above the melting line of the vortex lattice in the pristine 
crystal, H^T). This is at odds with the expectation that 
at = 90° the CDs act as weak point disorder; then, the 
melting and irreversibility lines are expected to lie below 
H^(T) 0. The angular dependence of H^ r (T) is more 
pronounced at lower temperatures. It decreases upon 
increasing T until it disappears completely at ~ 82 K, 
above which all lines coincide with H^(T). This behavior 
was observed for different samples with of 20, 30 and 
40 G. The relative shift in Hg r (T) between = to 90° 
was found to increase with B^. 

To examine the vortex behavior below H^ r (T) , within 
the vortex solid phase, we performed conventional MO 
measurements, in which a sequence of images is taken 



while the applied field H a is swept. Figure 1(c) shows 
flux profiles across the irradiated regions. These have a 
linear (Bean-like) shape due to the enhanced pinning in 
these regions. The profiles were taken at H z = 22 Oe and 
T = 67 K, well below H?? r (T). At = 90°, the profile is 
very shallow. In contrast, at 6 = 0, a clear sharp profile 
is observed. The critical current, estimated from the flux 
gradient, is 3.5 times higher at — than at 90°. 




FIG. 2: (color online) The melting transition in crystal B. (a) 
DMO image of part of the sample at H z = 35 Oe, T = 81 K. 
(b) Local DMO intensity versus T for different angles 0. The 
peak around 80.4 K corresponds to the melting transition. 
The curves are vertically displaced for clarity. 

The angular dependence of the melting transition was 
studied using DMO upon temperature sweeps with T 
modulation jl], [l8j]- The first-order melting transition 
is characterized by a step in the equilibrium flux den- 
sity, and appears as a bright signal in the DMO images. 
For example, the bright circular regions in Fig. 2(a) are 
the irradiated areas that undergo melting, while the sur- 
rounding pristine regions are already in the liquid phase 
and hence appear as a darker background. 

Figure 2(b) shows the local DMO intensity, averaged 
over a region of 15 x 15 /i m 2 in one of the irradiated 
regions (red square in Fig. 2(a)). The melting transi- 
tion manifests itself as a peak in the local intensity as T 
is swept. Each curve corresponds to a different 0, while 
H z = 35 Oe was kept constant. In contrast to the irre- 
versibility line that shifts to lower T upon increasing 0, 
the melting temperature remains unaffected, even when 
B is perpendicular to the CDs. Moreover, the transition 
width, determined from the width of the peak, appears to 
be the same for all measured angles. Similar results were 
obtained for different H z values. Note that the pristine 
melting line was found to have a weak linear dependence 
on the in-plane field [20] which cannot be detected within 
our experimental resolution for our low in-plane fields. 

Fig. 3(a) gathers the melting lines obtained at differ- 
ent angles 6. They all fall on a single line, H^(T), that 
terminates at a critical point (CP) below which the tran- 
sition becomes continuous in nature. The CP occurs at 
the same temperature of 77 K at all angles. The width of 
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FIG. 3: (color online) (a) The melting line H^{T) and the 
derealization line HdiiT) in crystal A, mapped for different 
0: 0(+), 25° (o), 45°(#), 65° ( A) and 90° (x). (b) DMO self- 
field image of part of sample A at 40 Oe, 75 K, 20 mA. 

the transition, denoted by the horizontal right-pointing 
arrows, is broadened on approaching the CP, but remains 
independent of angle, as was already shown in Fig. 2 for 
H z = 35 Oe. In Fig. 3, we also plot the derealization 
line HdiiT) Q, which separates the nanoliquid from the 
homogeneous liquid. HdiiT) was determined using MO 
visualization of the self-field induced by a transport cur- 
rent applied to the sample. DMO image of the self-field 
allows one to visualize the transport current density dis- 
tribution within the sample. Below HdiiT), the current 
flow is highly nonuniform due to enhanced pinning in the 
irradiated regions [2], and the DMO image shows signif- 
icant spatial variations (Fig. 3(b)). These disappear at 
Hdu when the network of pinned vortices is delocalized, 
and a homogeneous liquid is formed. As H^(T), the 
HdiiT) does not depend on the field angle (Fig. 3(a)). 
On approaching T c , vortex pinning by CDs becomes weak 
and hence all the transition lines coincide. 

The experimental observations show that the dynamic 
properties of porous vortex matter are strongly depen- 
dent on the angle between B and the CDs, while the ther- 
modynamic melting and derealization lines are angle- 
independent. In order to understand these observations, 
we have carried molecular dynamics simulations of a rect- 
angular parallelepiped model system of 36 PV stacks and 
200 layers. As shown in 0, [2l[, such system size is suf- 
ficient for observing the various phase transitions. Peri- 
odic boundary conditions (PBC) were implemented in all 
directions. The intraplane PV interaction was modeled 
by the repulsion of their screening currents, while the 
interplane PV interaction is described by their mutual 
electromagnetic attraction. For details see Refs. 0,|23. 
We also included the attractive Josephson interaction be- 
tween pairs of pancakes in adjacent planes, but belonging 



to the same stack [22|. We keep track of pancakes be- 
longing to a vortex stack and allow for flux cutting and 
recombination [22[. The average vortex direction (i.e. 
H a ) was kept parallel to the z axis, while tilted CDs of 
radius r r = 30 nm were introduced at random positions. 
The interaction between a pancake and a CD in the same 
plane is given by an attractive potential 0, [ll[ , which has 
a long range tail contribution « —e^dr^/R 2 , due mainly 
to the electromagnetic pinning, and a short range flat 
region of depth « e$d due to the combined effects of 
electromagnetic and core pinning. Here, eo = ((/>o/47rA) 2 , 
A = A(0)/(1 — T/T c ) 1//2 is the penetration depth, R is the 
distance between a pancake and a CD and d is the layer 
separation. The interaction — B • H/(47r) is implemented 
by the PBC, since the electromagnetic interaction be- 
tween a tilted stack and the stacks of images, positioned 
above (and below) its center of mass, pull the tilted stack 
in the z direction. Temperature effects are implemented 
in the simulation by using a white thermal noise whose 
variance is proportional to the temperature. 

Figure 4 shows projections of the vortices onto the xy 
plane (top view) and the xz plane (side view), with CDs 
tilted by 45° and 80°. The main observation is that, for 
B > B<f,, the CDs are almost fully decorated by pancakes 
both in the solid and nanoliquid phases. This means that 
the matrix of pinned PVs is preserved even when the field 
is tilted away from the CDs, and as a result the intersti- 
tial PVs remain caged in the surrounding fixed matrix. 
The interstitial PV stacks are even found to be tilted to 
some extent in the direction of CDs, in the solid phase. 
Another important observation is that the vortices adopt 
a structure featuring kinks between the layers, as shown 
by the side views in Figs. 4(c) and 4(d). Kinks of dif- 
ferent vortices tend to align, forming a structure which 
resembles the Josephson and Abrikosov crossing lattices 
(2^ in BSCCO at tilted field. Here, however, the Joseph- 
son vortices (JV) are formed when tilting the field away 
from CDs rather than from the c-axis. For example, when 
B is parallel to the CDs but at 45° with respect to c-axis, 
no JVs are present, while for B || c-axis JVs segments ap- 
pear, as shown in Figs. 4(c) and 4(d). 

To further interpret the results we consider the energies 
involved: The binding energy of a vortex to a CD is of 
the order of e^d per layer. For angles less than about 75° 
this gain exceeds the electromagnetic energy loss. For an 
infinitely long stack tilted at angle the later is given by 
eodln((l + cos(0))/2cos(0)) per layer [24]. In addition, 
the energy — B • H/(47r) tends to align the average vor- 
tex orientation with H. The contradicting requirements 
imposed by the need to minimize all these energies are 
optimized by the creation of kinks. The Josephson en- 
ergy cost of the kinks is negligible compared to the total 
energy of the vortex system in a highly anisotropic ma- 
terial like BSCCO jlHj. Indeed, for 7 - 375, the total 
energy measured in the simulations for non tilted CDs 
and for CDs tilted at 45° is the same, in the vicinity of 
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FIG. 4: (color online) Snapshots of vortex stacks and CDs 
tilted at 45° (a, b and c) and 80° (d) for B = 100 G and 
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35 G. Pancakes in the same stack are connected. Free 



pancakes are in blue, trapped ones in red. Only CDs at the 
bottom layer are shown (green). (a)(b) Projection onto the 
a-b plane (top view), (a) Nanosolid phase, (b) Nanoliquid 
phase, (c) Projection of (a) onto the a-c plane (side view), 
first row is shown, (d) Side view with CDs at 80°. 



the melting transition (T = 73 K to 77 K), within the 
simulation's error bars. The comparison was done for the 
same value of B e ^ = cos 6 in order to match the ex- 
perimental situation where the ratio of vortices to CDs 
remains fixed as 6 changes. Our simulations show that 
the rigid matrix of pancakes pinned on CDs is preserved 
up to very large angles in both the solid and nanoliquid 
vortex phases. This means that the enhanced caging po- 
tential created by the CDs, that was found to shift the 
melting line upwards, [1] is negligibly affected by increas- 
ing the angle between B and the CDs. Consequently, the 
thermodynamic properties, in particular the melting and 
derealization lines, are angle independent. This was ver- 
ified in our simulations for B = 100 G, where the melt- 
ing and derealization transitions occur at 74 K and 76 
K respectively, both for nontilted and tilted CDs, as ev- 
idenced by measurements of the mean square deviation 
of the flux lines, the amount of their entanglement and 
other measurements [26]. 

These observations can resolve the apparent discrep- 
ancy between the dynamic and thermodynamic proper- 
ties. The kinks developed by the vortices have a neg- 
ligible effect on the thermodynamics but a large effect 
on the dynamics. As demonstrated by our simulations, 
the kinks slide along the CDs due to the Lorentz force 
exerted by in-plane current, see also [l(J Conse- 
quently, they transfer PVs between segments, thus pro- 
ducing an effective transverse drift of the whole pancake 
stack. Upon increasing the amount of kinks increases 
and the flux creep becomes more effective. This process 
decreases the critical current and shifts the irreversibility 
line to lower temperatures. This is consistent with [141 ]. 



which concluded that CDs' pinning anisotropy commonly 
measured in BSCCO [17| must be a dynamic effect. 

In summary, the porous vortex matter in BSCCO is 
shown to preserve its thermodynamic properties even 
when the field is tilted away from the CDs. While the ir- 
reversibility line shifts to lower temperatures upon tilting 
the field, the thermodynamic melting and derealization 
lines remain constant. Numerical simulation shows that 
increasing the angle between the field and the CDs leads 
to formation of weakly pinned vortex kinks, while pre- 
serving the basic structure of a rigid matrix of pancake 
stacks residing along the CDs with nanocrystals of inter- 
stitial vortices embedded within the pores of the matrix. 
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